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ABSTRACT 

Emission line profiles of tracer molecule H 2 CO 140 GHz transition from gravitational 
core collapsing clouds in the dynamic process of forming protostars are calculated, 
using a simple ray-tracing radiative transfer model. Three self-similar dynamic inside- 
out core collapse models - the conventional polytropic model, the empirical hybrid 
model and the isothermal model - for star-forming molecular clouds are examined and 
compared. The isothermal model cannot produce observed asymmetric double-peak 
molecular line profiles. The conventional polytropic model, which gives flow veloc- 
ity, mass density and temperature profiles self-consistently, can produce asymmetric 
double-peak line profiles for a core collapsing cloud. In particular, the blue peak is 
stronger than the red peak, consistent with a broad class of molecular line profile 
observations. We find that line profiles are robust against variations in the polytropic 
index 7 once the effective line-centre opacity kq is specified. The relative strengths of 
the blue and red peaks within a molecular line profile are determined by the cloud 
temperature gradient, but the emission at frequencies between the two line peaks is 
determined by detailed density and velocity profiles in the cloud core. In the presence 
of a static dense kernel at the centre of a collapsing cloud, strong internal absorp- 
tion along the line-of-sight may occur, causing a suppression to the red wing of the 
blue line peak. If reliably resolved in frequency by observations, this signature may 
be potentially useful for probing the environs of an infant protostar. The conventional 
polytropic model can be utilized to produce molecular line-profile templates, for ex- 
tracting dynamical information from line spectra of molecular globules undergoing a 
gravitational core collapse. We show a sample fit using the 140 GHz H2CO emission 
line from the central region of the molecular globule B335 by our model with 7 = 1.2. 
The calculation of line profiles and fitting processes also offer a scenario to estimate 
the protostellar mass, the kernel mass accretion rate, and the evolution time scale of 
a core collapsing cloud. Our model can be readily adapted to other tracer molecules 
with more or less constant abundances in star-forming clouds. 
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— line: profiles — radiative transfer — stars: formation 



1 INTRODUCTION 

How molecular clouds collapse to form protostars is an 
important question in modern astronomy and cosmol- 
ogy. Since hydrodynamic collapse models for star forma- 
tions in_j]3rjle£uljj£^ proposed in the 1960s 
(e.g. iBodenheimer fc Swe igart 19681; iLarsonl [l969l : iPenstonl 
1 19691 ). substantial progress has been achieved in mod- 
elling dynamic properties of core collapsing clouds (e.g. 
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20091 : lLou fc ZhaJ 120091 ) . These prog: 



sive research advancement leads us to a more realistic un- 
derstanding of the star formation processes in molecular 
clouds. In particular, the gas flow velocity, mass density 
and thermal temperature structures of collapsing clouds 
predicted by theoretical models are in principle testable 
by various high-resolution spectral imaging observati ons in 
differ e nt bands of trace r molecular line profiles (e.g. Snelll 



1 19811 : IZhou et all 1 19951 : iMardones et~ail 1 19971 : lYang et all 
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20021; |Pe Vries Narayanan fc SnelH 120021 ; ISaito et all 1 19991 ; 
Belloche et al.ll2002l ; lEvans et alj|2005h . 



A widely discussed scenario for the formation of low- 
mass protostars is the so-called 'inside-out collapse' model 
((ShuJll977D. In this model, a condensed cloud core becomes 
unstable gravitationally and collapses rapidly towards the 
centre to form a protostellar core. The protostar is embed- 
ded in a dusty envelope, which is presumed to be static, 
although in some variations of the mode l the envelope may 
also involve an in f all or an outflow (e.g. Lou fc Shenl 12004 ; 
IShen fc Loul 12004 iHu fc Loull2008l ; lLou fc Zhaill2009f ). The 
boundary between the collapsing core and the envelope ex- 
pands, engulfing more and more mass into the central col- 
lapsing core. The accretion of dusty materials from the en- 
velope forms a disc rotating about the central protostar, 
and magnetized bipolar outflows may occur from th e disc 
at the same time (e.g. IShu. Adams fc Lizand [l987). The 
inside-out model has enjoyed supports from observations of 
some molecular globules. The B o k mo le cule globule B335 



(e.g. iFrerking. Langer fc Wilsonl [l987l; IZhou et al l 1 19931 : 



IChoi et al.lll995l ; ISaito et al.lll999l ; IStutz et al.ll200Sf) serves 
as a t est example among others (e.g. lHogerheiide fc Sandelll 
2000). However, there are features of li ne emissions 
from this star forming cloud (see e.g. IWilner et al.l 
2000) and also observati ons of other cloud system s, such 
as IRAM 04191 ( e.g. I Andre. Motte fc Bacmannl " 



Belloche etahl |2002|) and L1544 (e.g . iTafalla et all 



1999 



1998 



van der Tak. Caselli fc Ceccarelhf 2005), that are not readily 



explained by the isothermal inside-out collapse scenario as 
advanced by Shu et al. (1987) 

Different from the inside-out collapse scenario for a 
molecular cloud, there are alternative models to describe the 
formation of low-mass protostars within molecular clouds, 
e.g. the isothermal Bonnor-Ebert sph ere embedded in an 
exter nal medium of finite pressure P B l|Bonnorlll95l lEbertJ 
1955). This model describes the hydrostatic equilibrium of 
a gaseous sphere with a constant kinetic temperature, while 
a central core with finite density is involved. Such a hy- 
drostatic sphere can be maintained as long as the mass 
of the cloud does not exceed the critical value of M c = 
1.18a 4 /(G 3/2 Pe /2 )> where a is the isothermal sound speed 
and G is the gravitational constant. With density profile 
and kinetic temperature being the two key aspects, dynamic 
instability is predicted therein, but no flow velocity is pre- 
sumed. Recent research shows that the hydrostatic Bonnor- 
Ebert sphere may match several observed properties of star- 
less c ores (e.g. lEvans et alj|200ll ; [Belloche et aI1l2002l ; lMversl 
120051 ). while it cannot account for the observational data of 
B335 source (see e.g. lHarvev. Wilner fc Myers! 120031 ). This 
isothermal sphere model with a finite mass in a hydrostatic 
equilibrium cannot fully describe observations of molecular 
line profiles in the cloud core because of the very isother- 
mal assumption as we shall show presently in the following 
sections. 

T he Larson-Pen ston (LP) self-similar solution (|Larsonl 

1 19691 ; iPenstonl Il969h of an isothermal sphere is a clas- 
sic dynamic model for low-mass star formation. This LP 
type model predicts the formation of stellar density cores 
around the centre of molecular clouds. Characterisitcally, 
the decreasing-towards-centre radial collapse speed profiles 
is obtained along with corresponding mass density pro- 
files. However, the isothermal assumption encounters the 



similar problem of the Bonnor-Ebert sphere for modelling 
the observed molecular line profiles. Relaxing the dynamic 
self-consistency, all these models may be partially applied 
semi-empirically to star formation clouds for developing di- 
agnostics. For example, recent empirical works make use 
of the millimetre continuum emission and provide em- 
piric al thermodyna mic ( density and temperature) p rofiles 



_ Tafalla, et ail |2002| ; IShirlev. Evans fc Rawlingsl |2002| ; 
lEvans et al.ll2005f T These model results offer an independent 
perspective to probe star-forming clouds, but the overall dy- 
namic profiles lack the desired self-consistency. In reference 
to these dynamic models, current observation s and empiri- 
cal works, the polytropic EECC solution ( e.g. lLou fc Wand 
l200rj ; lLou fc Gadl2006l ; IWang fc LovfeOOSl , Lou & Hu 2009) 
appears to be a sensible alternative. With radial velocity 
fields, corresponding density profiles and variable tempera- 
ture profiles, this self-similar dynamic model should be more 
practical and versatile in explaining different observations 
(Lou & Gao 2006). 

Molecular emission lines provide valuable observational 
diagnostics to probe the dy namics and physical cond itions 
in star- forming clouds (e.g. iDvson fc Williams! 1 19971 ). For 
examples, optically thick emission lines of tracer molecules 
H 2 CO, HCO+, CO, CS and N 2 H+ are frequently identified 
for such diagnostic purposes. However, there are still un- 
known and uncertain parameters as well as model degenera- 
cies to explain the line strengths and profiles, due to global 
hydrodynamics, macro- and micro-scale kinematics, thermal 
structures, chemical compositions and turbulence. It is thus 
important to resolve these issues step by step to avoid ambi- 
guities in spectral modelling as well as in both observational 
and theoretical interpretations. 

To date, most molecular line profile calculations assume 
a certain parametrization for the density, velocity and tem- 
perature st r uctur es within emitting clouds. For example, in 
IZhou et al.ll| 19931 ). the so-called 'exact theory' takes the den- 
sit y and veloc ity profiles from the isothermal dynamic model 
of IShul (|l977f ) but the temperature profile is assigned sep- 
arately with a certain empirical base and therefore the gas 
cloud model is actually non-isothermal for molecular line 
calculations. This is a worthwhile attempt, yet such a speci- 
fication of the radial temperature structure is unsatisfactory 
as there are inconsistencies in the velocity and density pro- 
files of an isothermal cloud from the t heoretical perspec- 
tive. In recent years, calculations (e.g. |Pe Vries fc Myers] 
2005) considered more structural complexities in star- 
forming molecular clouds and invoked a more sophisti- 



evant atomic processes (e.t 


1. IWard- Thompson & Buckley 


2001 


Rawlings & Yates 2001 


; Keto et al. 2004; Tafalla et al. 


2006 


Tsamis et al.l2008; Pavlvuchenkov et alj|2008f). These 



studies aimed a t resolving specific issues, such as the chemi - 
cal profiles fe.g. lRawlings fc Yatedl200ll ; lTsamis et al.ll2008l ) 
or the competition betwee n relevant atomic processes (e.g. 
iPavlvuchenkov et al.ll2008l ). As the effects of dynamical and 
thermal structures were not their prime objective, most of 
these calculations basically parameterized the density, ve- 
locity and temperature structures in a collapsing cloud or 
simply adopted the isothermal inside-out collapse model of 
Shu (1977) but took an empirical temperature profile from 
observations of dust continuum. 

The key question now is: does it make a difference if the 
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structure profile of a cloud is determined self-consistently 
through a polytropic hydrodynamic model analysis instead 
of using a simple isothermal inside-out model but with a 
non-isothermal temperature profile? Also, are the molecu- 
lar line profiles different for models with different polytropic 
indices of the gas dynamics? These questions are essential 
for our understanding of the dynamics and physical condi- 
tions, in particular, in the early stage of forming low-mass 
protostars. The objective of our study here is to address 
this issue quantitatively. For instance, we calculate H2CO 
140 GHz molecular line profiles using a self-similar hydro- 
dynamic model that involves a conventio nal polytropic equa- 
tion of state (EoS) (Fatuzzo et al. 2 004; iLou fc Wanell2006l ; 
lLou fc Gaol200rj ; IWang fe Loul200Sl , Hu & Lou 2 008) and an 
isothermal inside-out collapse model (|Shulll977h . and care- 
fully compare their differences. We also directly compare 
the H2CO 140 GHz molecular line profiles obtained from 
the conventional polytropic model (|Lou fc Gad [2006) with 
observations, using the molecular cloud B335 as a test exam- 
ple. Comparison of our theoretical approach with empirical 
approaches used in former works is also made. In principle, 
our model calculations can be readily adapted to line pro- 
files of tracer molecules of more or less constant abundances 
in star-forming clouds. 

This paper is structured as follows. In Section 2, we 
present model structures of conventional polytropic and 
isothermal inside-out collapse models with their flow veloc- 
ity, density and temperature profiles. The relevant parame- 
ters and scalings of dynamic models are also specified there. 
We provide an overview for the formation of molecular line 
profiles in star-forming clouds in this section. In Section 3, 
we present the results of calculations, the comparison of the 
lines obtained by dynamical models and the comparison of 
the theoretical with the observed molecular line profiles. Ef- 
fects due to the cloud temperature gradient, optical depth, 
intrinsic line width, polytropic index and the presence of a 
relatively dense central kernel are also considered. Finally, 
results of our analysis are summarized in Section 4. For the 
convenience of reference, we summarize the basic polytropic 
hydrodynamics and the relevant asymptotic self-similar so- 
lutions in Appendix A, elaborate on how the blue and red 
peaks of a molecular line profile are shaped in an optically 
thick cloud in Appendix B, and lay out the radiative trans- 
fer formulation and our ray-tracing computation algorithm 
in Appendix C. 
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Figure 1. Parameter variations as functions of polytropic index 
7 = 2 — n for the expansion-wave collapse scenario with static 
SPS outer envelopes and dynamic free-fall cores as in Model I. 
The upper panel is the infall radius xq, the middle panel is the 
mass parameter A given by eq. (f2j) for a SPS solution d A 11 I t and 
the bottom panel is the point mass mo for x — > + . Consistently, 
7 = 1 is the isothermal case (Model II) with parameters xq = 1, 
A = 2 and mo = 0.975. All these parameters increase slightly 
and then fall as 7 increases from 1 to (4/3) ~ . The special case of 
7 = 4/3 (n = 2/3) is investigated by Lou & Cao (2008) and Cao 
& Lou (2009). 



Table 1. Parameters of self-similar polytropic dynamic solutions 
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n 




mo 


A 
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1.0 


1.0 


1.00 


0.975 


2.00 





1.1 


0.9 


1.17 


1.25 


2.04 





1.2 


0.8 


1.35 


1.15 


1.66 





1.3 


o.v 


1.28 


0.26 


0.40 






In all these three models, the molecular cloud consists of a 
free-fall collapsing inner core and a static outer envelope. 
A self-similar dynamic solution is adopted for the collaps- 
ing core in both Model I and II. For Model I, we further 
use three specialized solutions for comparison. In fact, the 
isothermal case (Model II) can be regarded as a very special 
case of the conventional polytropic model (Model I) with 
7 = n = 1. 



2 SELF-SIMILAR HYDRODYNAMIC MODELS 

Three self-similar hydrodynamic models for the velocity, 
density and temperature profiles are considered in our 
molecular line-profile calculations, namely Model I, a hy- 
drodynamic coll apse model with a conventional polytropic 
EoS for the gas (|Lou fc Gadl2006l . Hu & Lou 20 08); M odel 
II, an isothermal inside-out collapse m odel (|Shrj 1977|); a nd 
Model III, an empirical hybrid model (jZhou et all 1 19931) fl 



1 The hybrid model of Zhou et al. (1993) takes the mass density 
and flow velocity profiles from the isothermal model of Shu (1977) 
and uses an empirical variable temperature profile. This empirical 
temperature profile was inferred based on an isothermal density 
profile which scales as r~ 2 . 



2.1 Parameters of Hydrodynamic Model Solutions 

Details of Model I and relevant asymptotic solutions are 
contained in Appendix A for reference. Model I is charac- 
terized by 71 + 7 = 2 where n is a scaling index for polytropic 
self-similar hydrodynamics and 7 is the polytropic index (see 
Appendix A; Lou & Wang 2006; Lou & Gao 2006; Hu & Lou 
2008). We choose different dynamic solutions with a static 
outer envelope and free-fall inner core collapse for compari- 
son in terms of tracer molecule line profiles. The density, ve- 
locity and temperature profiles can be directly derived from 
equations (|A9|I . (IA10|) . (|A7|) for a conventional polytropic 
gas dynamics under self-gravity. 

For conventional polytropic solutions with 7 = 1.1, 1.2 
and 1.3, the corresponding values of the dimensionless infall 
radius xo, the dimensionless core mass mo, and coefficients 
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Figure 2. Inward radial flow velocity, number density and ther- 
mal temperature profiles (panels from top to bottom respectively) 
of three different model molecular clouds used in our molecular 
line-profile calculations. The abscissa represents the radius R in 
unit of 10 3 AU, shown on a logarithmic scale. The solid curves 
denote Model I with 7 = 1.2; the dashed curves denote Model 
II. The dash-dotted curve in the bottom panel denote Model III, 
whose redial velocity and number density has the same manner 
as Model II (dashed curves in the top and middle panels) . Other 
relevant model parameters can be found in Table 1 of 7 = 1.0 for 
Model II and III, and 7 = 1.2 for Model I. The dimensional 
scalings used are described in Section 2, giving infall radii of 
R = 10 X 10 3 AU for Model II and III, and R = 13.5 X 10 3 
AU for Model I (7 = 1.2), respectively. 



A and B for asymptotic solution (|A12|l are contained in Ta- 
ble Q] Parameters for the isothermal model of 7 = n = 1 are 
also listed in Table[T] In our constructed solutions with static 
singular polytropic solution (SPS) envelopes, the dimension- 
less infall radius xo depends on n only and is analytically 
given by 



x 



(2- 



Nl/2 



[2(3n-2)] 1 "" /2 



(1) 



The velocity parameter B of asymptotic solution for static 
envelope vanishes and the mass parameter A can be derived 



A = 



2(3n - 2) (2 - n) 



1/r. 



(2) 



which is exactly the same as the coefficient of a(x) in SPS 
solution (|A11|) . The dimensionless central point mass mo of 
asymptotic free-fall solution (IA13[) is then determined nu- 
merically. Fig. [T] shows parameter variations of xo, A and 
mo as the scaling index 7 (or equivalently n = 2 — 7) varies 
continuously. This is in fact the case (i.e. expansion-wave col- 
lapse solution: EWCS) of a static SPS envelope connected 
with a free-fall collapsing core an d is much more general 
than the special isothermal case of Ishul (|1977| ). 

The dimensionless variables are converted to corre- 
sponding dimensional physical variables in order to perform 
radiative transfer calculations for molecular line profiles in 
core collapsing clouds. We adopt parameter values estimated 
from observations for the scalings of the dynamic profiles for 
model clouds. The typical infall radius within a mole cular 
cloud is estimated as ~ 0.01 — 0.03 pc (e.g. [Myers 2005j), cor- 
responding to ~ 10 3 — 10 4 AU. As dimensionless infall radii 



are roughly around and larger than xo = 1 in our model 
(e.g. 1 < xo < 1.4 by Fig. [T]and Table 1), we may choose 
the following length scale 

k 1/2 t n = 10 4 AU (3) 

in definition ()A5|) . A typic al number density at the infall 
radius is ~ 10 4 c m" 3 (e.g. lHarvev. Wilner fc Myersl 120031 ; 
lEvans et al.ll2009l ). For the mass parameter A < 2 as shown 
in Fig. [T] the reduced number density being of order of unity 
implies a proper number density scale as (see equation IA6[) 



l/(47rG/M7i H r) = 10 4 cm" 



(4) 



This corresponds to a thermal temperature of T ~ 4 K by 
equation (|A7|I and a characteristic infall age of t ~ 3 x 
10 5 years b y equation Q for a molecular cloud (see e.g. 
lMversll2005h . Within the infall radius, this scaling also gives 
an enclosed mass scale of M ~ 0.5 M© and a sound speed 
scale c ~ 0.15 km s^ 1 by equations ()A6|I and (IA8|) . 

In addition to the central protostar, we presume a num- 
ber density of ~ 10 6 cm -3 at the very central region (i.e. 
R < 500 AU) of a molecular cloud and set the central 
velocity to be zero for a static dense kernel. The outer 
cloud radius is set to be several tens of 10 3 AU. We note 
that the size of the static dense kernel and the outer ra- 
dius of the clo ud are not the same as used in some previ- 
ous works (e.g. lHarvev. Wilner fc Myersl 120031 ; iMverdbooj ; 
lEvans et al.ll2005r i. The key parameter in our radiative trans- 
fer calculations is the effective optical depth of a molecu- 
lar line in a cloud. Once the numerical value of the optical 
depth along the ray is specified, the actual size of an emit- 
ting cloud would be irrelevant. In this sense, the difference 
between cloud sizes used in our model analysis and those 
of earlier research works should not affect our basic conclu- 
sions. Nevertheless, we should keep in mind that this con- 
clusion is made under the assumption that the abundance 
of tracer molecule remains constant in our model consider- 
ation. The cloud sizes would matter if the fractional abun- 
dances of certain tracer molecules vary considerably in star 
forming clouds. 

For comparison, the velocity, number density and tem- 
perature profiles of the three different models that we use in 
radiative transfer calculations of molecular line profiles are 
shown in Fig. [5] In all models, the velocity and temperature 
increase monotonically towards the centre until R « 500 
AU. There are slight differences between the density and 
velocity profiles of Model I (7 = 1.2) and Model II. The in- 
fall radii for Model II and Model I (7 = 1.2) are at R = 10 
and 13.5 AU respectively. While the temperature profiles of 
the three models are readily distinguishable. Model II has a 
flat constant temperature, Model I (7 = 1.2) has a mono- 
tonic increase of temperature towards the centre, and Model 
III has a semi-empirical profile in which the temperature 
drops slightly from R = 30 x 10 3 AU to R « 6 x 10 3 AU 
and then increases when R decreases further. We also show 
in Fig. [3] the velocity, density and temperature profiles for 
Model I with different polytropic indices: 7 = 1.1, 1.2 and 
1.3, whose infall radii are at R = 11.7, 13.5 and 12.8 xlO 3 
AU respectively. As shown, cloud structures can vary sub- 
stantially even within Model I (in fact, Model II can be 
regarded as a special case of Model I), and such differences 
can leave signatures in the emission molecular line profiles 
as expected. 
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Figure 3. Inward radial flow velocity, number density and ther- 
mal temperature profiles (panels from top to bottom respectively) 
of Model I with three different polytropic indices 7. The dashed, 
solid, and dash-dotted curves represent 7 = 1.1, 1.2, and 1.3 with 
infall radii at R = 11.7 X 10 3 , 13.5 X 10 3 , and 12.8 X 10 3 AU, re- 
spectively. The corresponding central mass points are 1.13, 1.02 
and 0.24 Mq. Other relevant polytropic model parameters are 
contained in Table 1 and the physical scalings are described in 
Section 2. The case of 7 = 1.2 is the same as shown in Fig. [2] A 
larger 7 provides steeper density and temperature increases to- 
wards the centre; this effect will manifest in the corresponding 
H2CO 140GHz line profiles as discussed in Section 3.3. 



Figure 4. Profiles of H2CO 140 GHz emission lines in terms 
of brightness temperature Tb from molecular clouds for Model 
I with 7 = 1.2 (solid curves), Model II (dashed curves) and 
Model III (dash-dotted curves). The radial flow velocity, parti- 
cle number density and thermal temperature profiles of the un- 
derlying dynamic models are shown in Fig. [2] with a central ker- 
nel to avoid the singularity. The intrinsic line widths A» are set 
to be ~ 0.3 km s _1 here. Panels from left to right correspond 
to line emission profiles along the light-of-sight passing through 
R = AU, R = 1.2 x 10 3 AU and R = 6.2 x 10 3 AU, respectively. 
The effective total line-centre opacity is kq = 7 X 10 3 cm 2 g _1 in 
all three panels. 



2.2 A Molecular Line Profile Calculation 

We utilize a simple ray-tracing radiative transfer model 
(which does not involve an explicit treatment of atom and 
molecular transitions and their couplings with the radiation 
field; see Appendix C for details) to calculate the spectral 
profile of molecule H2CO (2i2-ln) emission line at frequency 
140 GHz. As the range of temperature variation in a cloud 
region is typically small (e.g. ~ 10 — 20 K), the variation 
in tracer molecular level po pulations is small a nd thus will 
affect line profiles little (e.g. lTsamis et al"1 l2008). 

It has been argued that gradients in distribu- 
tions of chemical elements as well as chemical evolu- 
tion are important in determining the strengths and 
profiles of molecular emission line s from star-formin 
clouds (e.g. iRawlings fc Yatesl l200ll ; iTsamis et ail I200" 



IPavlvuchenkov et al.l [2008|) . As we focus on the effect of 
dynamical cloud structures on tracer molecule line pro- 
files, the H2CO abundance is taken to be approximately 
constant in our model calculations. This simplification im- 
plies that our approach may be extended to those tracer 
molecules with more or less constant abundances in star 
forming clouds. This assumption of co nstant abundance 
was invoked in some earlier work s (e.g. IZhou et~al ] 1 19931 : 
IWard- Thompson fc Buckleyl l200ll ). while in others, the 
H2CO abundance in B335 cloud was assumed to be a step 
funct ion with a variation of ~ 10 times (e.g. lEvans et al.l 
2005). In fact, recent empirical works based on observations 
do show that several diagnostic tracer molecules (e.g. HO, 
H2CO, CN, CS) keep their abu ndance variations le ss than 
a factor of ~ 10 (e.g. Fig. 3 in ITsamis et alj|2008f ). Both 
theoretically and observationally, clarifications of the rela- 
tive importance of dynamical and chemical effects are essen- 
tial in using molecular line profiles as a reliable diagnostic 



probe to the core collapse dynamics in star-forming clouds. 
A more comprehensive description of star-forming clouds 
should involve self-consistent dynamic profiles and chemi- 
cal abundance profiles that are independently inferred from 
observations. 

In order to have an intuitive feeling for our model ap- 
proach, we demonstrate in Appendix B how the blue and 
red peaks of emission molecular line profiles are formed in 
an optically thick cloud. The radiative transfer analysis and 
the computing algorithm of our line profile calculations are 
described in Appendix C. The physical variables (i.e. tem- 
perature T, mass density p, projected velocity component 
v p etc.) of a core collapsing molecular cloud are determined 
from our polytropic hydrodynamic model (see Appendix 
A). The line-centre opacity Ko and the line width Av are 
set to values typical of spectral line observations of tracer 
molecules. In equations f|C 1 f) — (|C4|) . the mass density p is 
mainly that of the interstellar medium in a cloud. As a 
consequence, the ratio of molecular abundance is merged 
into the parameter Ko, which is then effectively a product of 
molecular line centre opacity and its mass fractional ratio. 

A simple assumption of local thermodynamic equilib- 
rium (LTE) is made for the radiative transfer calculation; 
this calculation does not account for scattering dominated 
processes, masers, or other non-LTE cases. 



3 RESULTS AND ANALYSES 

3.1 Temperature Variations in Molecular Clouds 

It is clear that Model II does not produce asymmetric 
double-peak H2CO 140 GHz line profiles (see the dashed 
curves in Fig. [4] and Equations (|B3j and (|B4[I in Appendix 
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B), an apparent contradiction to extensive line profile obser- 
vations in molecular globules. Therefore, unless compromis- 
ing the self-consistency in the thermal structure and dynam- 
ics to insert a semi-empirical temperature profile to replace 
the isothermal model, as in Model III, one cannot have a sen- 
sible comparison between the observed asymmetric double- 
peak line profiles of star-forming clouds and the line profiles 
calculated from Model II. In fact, equation ()B4[1 in Appendix 
B clearly demands that the temperature should increase 
towards the cloud centre in order to produce a molecular 
line profile with a blue peak being stronger than the red 
pe ak for a core col lapsing cloud. Fig. U (see also Fig. 9a 
of lZhou et ailll993l ) demonstrates that asymmetric molec- 
ular line profiles with a stronger blue peak emerge when 
Model III and Model I (7 = 1.2) are actually introduced. 
Recent radiative transfer works use empirical temperature 
profiles derived from the sub-millimeter continuum emission 
intensity are able to produce the observed asymmetric line 
profi l es (e.g. Ward-Thompson fc Bucklevll200ll ; lEvans et al.l 
l200ll ; iTsamis et al. 20081 ). However, self-consistency of the 
dynamic profiles should be demanded for a more satisfying 
model construction. 

We emphasize that the temperature profile is deter- 
mined simultaneously once we derive mass density and flow 
velocity profiles from conventional polytropic hydrodynarnic 
equations. Instead, when one uses an empirical temperature 
profile for molecular line profile calculations, this empiri- 
cal temperature profile and the density and velocity profiles 
are generally not consistent with each other in the hydro- 
dynamic and thermodynamic sense. In other words, a gas 
cloud with dynamic profiles having no self-consistency does 
not exist according to hydrodynamics and thermodynam- 
ics. It should be also noted that the inference of empirical 
temperature profiles must also rely on several assumptions. 
Our goal is to develop theoretical models grossly satisfying 
available observational constraints. 

Our calculations show that in both Model I (7 = 1.2) 
and Model III, the asymmetry in line profiles can be more 
pronounced in the off-centre region in a cloud, provided the 
ray of emission along the line-of-sight passes through the 
collapsing core of the cloud (compare the left and middle 
panels of Fig. 3]). The relative strengths of the blue and red 
peaks in a molecular line profile show a stronger contrast in 
Model III than in Model I (7 = 1.2). This sharper contrast 
in two peaks from Model III is due to the fact that there 
are more variations in the specified semi-empirical temper- 
ature profile than those in the temperature profile of Model 
I (7 = 1.2) that we adopt. Model II produces symmetric 
single-peak molecular line profiles in the resolved spectra 
by our numerical explorations (see three dashed curves of 
Fig. ©, while Model I (7 = 1.2) and Model III could pro- 
duce single peak line profiles when the rays are along the 
line-of-sight passing through only the static outer envelope 
or when velocity variations are insignificant in the emission 
region as compared with the intrinsic velocity broadening 
Av of H 2 CO 140 GHz line. 



3.2 Line Opacity and Intrinsic Velocity Width 

The left panel of Fig. [S] shows how the two peaks in H2CO 
140 GHz line profile vary with the value of the line cen- 
tre opacity kq. With density distribution being fixed, the 
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Figure 5. Profiles of H2CO 140 GHz emission lines in terms of 
brightness temperature Tb from molecular clouds at R = 1.2 X 
10 3 AU for Model I of 7 = 1.2 with different values of the effective 
line-centre opacity k = 0.6 X 10 3 , 1.6 X 10 3 , 2.6 X 10 4 cm 2 g _1 
(left panel), and for different values of intrinsic line widths Av 
(0.1, 0.3 and 0.6 km s _1 ) with a fixed k = 1.6 X 10 3 cm 2 g" 1 
(right panel). Unless otherwise stated, the parameters adopted 
are the same as those in Model I in Fig. [2] 



opacity fto characterizes the optical depth along the ray 
parallel to the line-of-sight. For forming asymmetric blue 
and red peaks, molecular line profiles need to be optically 
thick, in addition to the presence of variable temperature 
and velocity structures within molecular clouds involving 
gravitational core collapse. For molecular lines that are op- 
tically thin, the observed emissions sample the entire col- 
lapsing core. A cloud structure of spherical symmetry for 
the near and far hemispheres inevitably gives rise to sym- 
metric molecular line profiles, in spite of the temperature, 
density and velocity variations in the emission region. More- 
over, as the line is optically thin, the line intensity is well 
below the source Planck function. Hence the line intensity 
becomes weaker, i.e. the corresponding brightness tempera- 
ture is lower than the thermal temperature of line emission 
regions (see equation IB4I and discussions thereafter for an 
intuitive explanation for the formation of a double-peak line 
profile). 

For comparisons of dynamic models, we have adopted 
that the intrinsic line width Av — 0.3 km s~ . There are un- 
certainties in the values of Av and we have chosen a typical 
value here. It is generally believed that the intrinsic width 
of a molecular line is mainly caused by thermal broadening 
and micro-turbulence. The thermal broadening is estimated 
to be at a level of ~ 0.1 km s _1 for a typical molecular 
cloud with a thermal temperature of ~ 10 K. The turbu- 
lent broadening is fairly uncertain. Theoretical estimates 
depend on the assumed driving mechanism s (e.g. iMacLowl 
Il999l ; iLi fc Nakamural 120061; lMatznerll2007l ). Some observa- 
tional stmchN3s(e^|Lars^^ that the probable line 
width Av ranges from a few ~ 0.1 km s _1 to about several 
km s . 

Typically, infall velocities v in cloud cores would range 
from ~ 0.1 km s" 1 to a few km s _1 , as adopted in the dy- 
namic models explored here. The characteristic double-peak 
signature of molecular line profiles from a collapsing cloud 
may be erased or smeared if the line suffers a severe turbu- 
lent broadening (i.e. large At;). Right panel of Fig[5]demon- 
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Figure 6. Profiles of H2CO 140 GHz emission lines in terms of 
brightness temperature Tb from molecular clouds for the three 
polytropic models shown in Fig. [3] The dash-dotted curve corre- 
sponds to the model with 7 = 1.1; the solid curve to the model 
with 7 = 1.2; and the dashed curve to the model with 7 = 1.3. 
The scale free infall radii are at xo = 1.18, 1.35 and 1.27 respec- 
tively. For all three panels, kq = 7 X 10 3 cm 2 g" 1 . 

strates that for Model I (7 = 1.2) with a maximum infall 
velocity of ~ 2 km s , we can produce a double-peak line 
profile for Aw < 0.3 km s _1 but not for Aw > 0.6 km s _1 . 
Extensive observations of double-peak line profiles in many 
molecular globules imply that the turbulence in their line 
formation region should not give rise to velocity broadenings 
significantly larger than ~ 0.3 km s _1 with a maximum in- 
fall velocity ~ 2 km s _1 . By a further numerical exploration, 
we find that the maximum allowed velocity broadening bears 
a linear relationship with the maximum infall speed; the val- 
ues shown here are typical for star-forming clouds. 

3.3 Polytropic Index 7 

One appealing aspect of Model I is that we can obtain 
self-consistent velocity, density and temperature profiles for 
molecular clouds within the overall inside-out collapse sce- 
nario. The essence of this model is the conventional poly- 
tropic equation of state characterized by a polytropic index 
7. The special isothermal case (Model II) just corresponds 
to 7 = n — 1. Now our question is: how robust is the molec- 
ular line profile morphology with respect to the variation of 
7? 

We calculate H2CO 140 GHz emission line profiles us- 
ing Model I with three different polytropic indices 7 = 1.1, 
1.2 and 1.3. These 7 values are representative of astrophys- 
ical situations between limiting regimes, e.g. the isothermal 
process of 7 = 1, and the adiabatic processes of 7 = 7/5 
and 5/3 for diatomic and monoatomic gases respectively. We 
found that the overall morphology of molecular line profiles 
is quite robust with respect to variations of 7, provided that 
the effective line-centre opacities in the calculations take on 
the same values (see Fig. |SJ. Note that n + 7 = 2 for a 
conventional polytropic gas, i.e. a variation of 7 leads to a 
variation of n, which is related to the scalings of hydrody- 
namic variables. There are however some small and subtle 
differences among predicted molecular line profiles. For in- 
stance, larger 7 tends to give more prominent peaks in a 



Figure 7. Comparison of H2CO 140 GHz emission line profiles 
in terms of the brightness temperature Tb along the light-of-sight 
passing through the cloud centre for Model I (7 = 1.2) (left 
panel) and Model III (right panel). The difference between the 
two model curves within each panel is that the solid curve is the 
same as the cloud as in Fig.[2]and dashed curve has the innermost 
zero-velocity condensed kernel removed. Other parameters of the 
'kernel-less' model clouds are the same as those in Model I and 
Model III in Fig. [2] In both panels, k = 7 X 10 3 cm 2 g _1 . 

line profile, which can be attributed to a much steeper gra- 
dient in the temperature increase towards the cloud centre 
(Fig. [3j). Also, for a larger 7 value, the molecular line peaks 
persist to larger projected radii away from the cloud centre 
(see the right panel of Fig. [6} . Nevertheless, these differences 
might not be readily distinguishable in observations, given 
the beam smearing in the spectral imaging observations and 
the presence of various noises in molecular line spectra. 

3.4 The Effect of a High-Density Central Kernel 

In theoretical models of gravitational core collapsing clouds, 
a static central dense kernel is introduced to avoid the sin- 
gularity. Phenomenologically, the radius of the kernel marks 
the interface of an 'embryo' stellar kernel with the accreting 
gas in a gravitational free fall. This high-density kernel is 
more a theoretical hypothesis as w as used in former works 
(e.g. IChoi et al.lll995l ; iMverj 120051 ) but it doe s fit the ob- 
serve d sub-millimeter intensity profile well (e.g. lEvans et al.l 
120011 ). 

The static high-density kernel could leave indelible sig- 
natures in molecular line profiles from the central core of 
a collapsing cloud. As the effective kernel opacity is large, 
molecular emissions from the far side (i.e. behind the dense 
kernel) would suffer strong absorptions at the rest-frame 
line-centre frequencies, depending on intrinsic line broaden- 
ing Av. In contrast, molecular emissions from the near side 
(i.e. in front of the dense kernel) is unaffected by the kernel. 
This naturally leads to the suppression of the red wing of the 
blue line peak. In the left panel of Fig.|4j molecular line pro- 
files manifest a dip feature in the red wing of the blue peak 
when the line-of-sight passes through the dense kernel. For 
rays along the lines-of-sight that avoid this central dense ker- 
nel, this dip feature would disappear (see the middle panel 
of Fig.|?]as an example). Fig.[7]compares molecular line pro- 
files from both Models I (7 = 1.2) and III with and without 
a central kernel. The difference in molecular line profiles for 
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the case with and without the dense static kernel are not 
negligible. Moreover, the presence of a dip makes the red 
peak appearing broader than the blue peak. However, this 
difference in line profiles is similar to the difference caused 
by different values of 7 in Model 1 (the left panel of Fig[6]), 
as well as the difference between Model 1 and Model III (left 
panel of Fig[4j| . So the actual cause of this phenomenon on 
profiles can only be solved after testing the temperature and 
density profiles of the cloud by some other means (e.g. the 
millimetre continuum). 

Potentially, such a line feature may be utilized to probe 
the 'embryo' kernel surrounding a protostar in its develop- 
ment and evolution [e.g. pre-protostellar cores (PPCs) with 
detectable dust continuum emission at millimetre and sub- 
millimetre wavelengths] before the eventual emergence as a 
'fully functional' star of nuclear burning. However in reality, 
whether or not such dips may be unambiguously resolved 
in the observed line profiles from molecular globules is de- 
termined by a few factors. In particular, it depends on the 
spectral resolution actually achieved by observations. 

3.5 Comparisons with Observed Line Profiles 

It is straightforward that self-similar solutions of Model I 
combined with radiative transfer calculations (Appendix C) 
can be applied to fit observed H2CO 140 GHz line profiles. 
Empirical estimates of physical variables at certain epochs 
are introduced for calibration and profiles of variables are 
derived from Model I self-consistently. As an example of il- 
lustration, we specifically model the 140 GHz H2CO emis- 
sion line from the Bok m olecular globule B3 35. We take the 
observed line profile from IZhou et al.l |l993) , where the line 
emission was extracted along the line-of-sight path towards 
the cloud centre. The collapsing core (i.e. the inf alling region 
aroun d the centre) radius is ~ 30" according to lZhou et al.l 
l| 1993ft . By adopting a n estimated distance of ~ 250 pc to 
the Bok globule B335 (|Tomita. Saito fc Ohtam1ll979l ). this 
gives a radius of ~ 7.5 x 10 3 AU for the infalling region. 
However, the uncertainty in this distance to B335 is fairly 
large (see Table 2 of Tomita et al. 1979) and this radius 
may vary by more than a factor of two. We take this radius 
to be ~ 1 x 10 4 AU for the collapsing core in our model 
calculations for B335 cloud. 

The same 7 = 1.2 solution of Model I as in Fig.[2]is used 
to model the structure of molecular line formation region. 
The maximum infall velocity is taken to be ~ 0.50 km s _1 at 
0.6 x 10 3 AU and the radial receding displacement velocity 



- 8.30 km s" 1 
0.07 km s~\ 



The intrinsic 
and the effec- 



ted of the entire B335 cloud is 
width of the molecular line is ' 
tive line-centre opacity is ko = 10 4 cm 2 g _1 . We find that a 
dense static inner kernel is needed in order to produce line 
profiles similar to the observed profiles. The kernel radius is 
~ 0.5 x 10 3 AU and the mean particle number density in 
the kernel is ~ 10 6 cm -3 . Fig. [5] compares the observed line 
profile and the line profile produced by Model I (7 = 1.2). 
The good agreement between our model line profile and the 
observed line profile clearly shows the viability of the conven- 
tional polytropic hydrodynamic model and the applicability 
of our line-profile radiative transfer calculations. 

It is possible to derive the central point mass Mo and 
the total enclosed mass contained within the central kernel 
region Mk in Model I used in the above data fitting. Using 




Figure 8. Comparison of the 140 GHz H 2 CO (2i 2 - In) line 
profile towards the central region of the Bok molecular glob- 
ule B335 with the theoretical molecular line profile obtained by 
Model I with 7 = 1.2 (dashed line). This abundance ratio of 
H2CO molecule is taken to be 7 X 10 — 9 and kq = 10 4 cm 2 g — 1 . 
The radial receding displacement velocity of the entire cloud is 
fd ~ 8.30 km s — 1 . The infall speed is ~ 0.27 km s — 1 at a ra- 
dial distance of ~ 10 3 AU away from the cloud centre. The gas 
temperature of the central kernel is ~ 3.6 K. 



the mo value (see Table [T] and Fig. [TJ and scaling relations 
([3]) and in our dynamic model with self-similar trans- 
formation (IA6j) . we obtain the value of central point mass 
Mo = 1.O25M0, referred to as an estimate for the mass of 
th e central protostar ; this estimate is about t wice the value 
of IZhou et all l| 1993ft and IStutz et alJ (|2008f l a nd is ~ 10 
times larger than the mass value estimated by ISaito et al.l 
(1999). When we use the reduced enclosed mass m(x) at 
the boundary of the central kernel, the dimensional enclosed 
mass contained within the kernel region is derived to be 
Mu = 1.027Mq. Estimates of all these masses are not sen- 
sitive to the kernel radius but are mainly affected by the 
variation in the infall radius. Therefore, all these masses 
would have an uncertainty of a few times. It is also possible 
to derive the mass accretion rate in the periphery around 
the kernel, M h = 0.8 x 1O~ 6 M yr _1 . In this case, the 
timescale of a protostar formation can then be estimated 
as t = M k /M h ~ 1.25 x 10 6 yr, whic h is ~ 10 times t he 
infall timescale estimated in Model III (|Zhou et al.lll"993T ). 



4 CONCLUSIONS 

We have investigated and analyzed effects of radial dynam- 
ical and thermal structures of an inside-out core collapsing 
cloud with spherical symmetry on forming molecular line 
profiles emitted from star-forming molecular clouds. A sim- 
ple radiative transfer model employing a direct ray-tracing 
algorithm along the line of sight is constructed for computing 
H2CO 140 GHz emission line profiles in order to compare 
with those actually observed in a star-forming cloud. We 
model the gravitational core collapse of a molecular cloud by 
the self-similar hydrodynamics of a conventional polytropic 
gas under self-gravity in a self-consistent manner (Model I). 
We emphasize that an isothermal inside-out collapse model 
(Model II) cannot produce asymmetric molecular line pro- 
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files. In contrast, a variable radial temperature structure, 
which is crucial for determining asymmetric molecular line 
profiles, can be derived self-consistently together with the ra- 
dial velocity and density structures in Model I, which is dif- 
ferent from the semi-empirical approach (Model III) . We ex- 
plicitly demonstrate that Model I can indeed produce asym- 
metric double-peaked line profiles with stronger blue peaks 
than red peaks as observed in many star-forming molecular 
globules. We assess the relative importance of effects due to 
line-centre opacities, intrinsic line widths and polytropic in- 
dex 7 in Model I. Our radiative transfer calculations further 
show that the presence of a dense static inner kernel around 
the centre of a core collapsing cloud gives a characteristic 
signature of a suppressed red wing of the blue line peak in 
the molecular line profile. We show an application to line 
profile observations by modelling the 140 GHz H2CO emis- 
sion line profile from the centre of the Bok molecular globule 
B335. We note that sensible fits to the H2CO emission line 
profile of B335 requires a dense static kernel in addition to 
conditions of a polytropic core collapse model and an esti- 
mated line-of-sight receding speed of the entire molecular 
cloud. Using Model I and the molecular line profile fits, we 
could estimate the protostellar mass, the mass inside the 
static condensed kernel, the mass accretion rate just out- 
side the kernel radius, and the time scale of the protostar 
formation. 
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APPENDIX A: HYDRODYNAMICS OF A 
CONVENTIONAL POLYTROPIC GAS 

For the hydrodynamics of a conventional polytropic gas 
model of spherical symmetry without diffusive effects, we 
use the basic nonlinear fluid equations in spherical polar co- 
ordinates (r, 6, (j>), namely 



dp . 1 d . 2 \ „ 
-di + -^ r PU) = ° 
du du _ 1 dp 
dt dr p dr 

dM_ +u dM_ _ Q 
dt dr 



GM 



dM 
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Here the mass density p, radial bulk flow velocity u, ther- 
mal pressure p, and enclosed mass M are all functions of 
radius r and time t, and G = 6.67 x 10 -8 g _1 cm 3 s~ 2 is 
the gravitational constant. Equations (|A1|I and (|A2|) are the 
mass and radial momentum conservations; equation ()A3|) 
represents another form of mass conservation; the Poisson 
equation relating p and the gravitational potential <E> is con- 
sistently satisfied by the above equations. 
The polytropic equation of state reads 



p = K{t)p~< 



(A4) 



where 7 is the polytropic index and K(t) is a coefficient. 
The self-similar dynamics involves an independent variable 

x = r/{k 1/2 t n ) , (A5) 

where n is a scaling index and k is referred to as the sound 
parameter, and the corresponding self-similar transforma- 
tion of other relevant hydrodynamic variables are 

k 3 ' 2 t An - 2 m{x) 



a{x) 
(47rGt 2 ) ' 

a = k 1/2 t n ~ 1 v(x) 



M = 



(3n - 2)G ' 
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where a(x), m(x) and v(x) are dimensionless reduced mass 
density, enclosed mass and radial flow velocity, respectively. 
In reference to equation of state (|A4|) . we identify K(t) = 
fc(47rG) 7 ~ 1 t 2(7+n_2) . For n + 7 = 2, K becomes a constant 
for a conventional polytropic gas (e.g. Suto fc Silk Il988l ; 
lLou fc Wan3l2006l ; lLou fc Gaoll2006l ; iHu fc Loull200a ). The 



gas temperature T and the local polytropic sound speed c 
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where fee, p. and ran represent Boltzmann's constant, mean 
molecular weight and the hydrogen mass, respectively. For 
typical star-forming clouds, we shall take p = 1 for sim- 
plicity. This self-similar transformation deals with an im- 
portant subset of nonlinear solutions to polytropic hydro- 
dynamic partial differential equations (PDEs) (|A1[I — (|A4[) 
and provides the dynamic structures for radiative transfer 
calculations of molecular line profiles. 

With self-similar transformation ()A6|) in nonlinear 
PDEs (|A1|I — (|A4|) . we obtain two coupled nonlinear ordinary 
differential equations (ODEs) 



[(*• 



■ 7a 



7-1 



dv (v — nx) 2 



dx (3n - 2) 



2(v — 2)70: 7 1 . . 

H — — (v - nx)(n - l)v , 



x 



r. .2 7-11 da 

(v — nx) — 7a — — 



(«■ 



da 
adx 
nx)a 



= (n — l)v 
2(v — x)(v — nx) 



(3n - 2) 



(A9) 



(A10) 



By setting v — for all x > 0, we obtain a static singular 
polytropic solution (SPS) from equations (|A9|) and (|A10[) 







[2 7 (3n - 2) jn 



21 V" -2/n 
X 



(2 7 ) 1/n (3n-2) 1/n n- 7/n a; (4 - 37)/ " . (All) 



To leading orders, asymptotic similarity solutions of coupled 
nonlinear ODEs JA§J and (|ATo)l for x — ► +00 and x — > + 
are summarized below. In the limit of x — > +00, we have 



a = Ax 



-2/r. 



nA 



(3n-2)' 



„(n-2)/n 



+ 



27A 



7-1 



(2-2 7 -n)/n 



+ Bx 



(n-l)/n 



(A12) 



where A and B are two constants of integration. In asymp- 
totic solution (IA12|) . it is possible to set B = with a finite 
A > such that v vanishes at large x; and this A is precisely 
the coefficient of a(x) in SPS solution (|A11[) . This consis- 
tency is expected but SPS solution (|AllJl is valid for all x 
except for the central singularity. In the limit of a; — + + , 
the asymptotic central free-fall solution is 



2mo 



(3n - 2)x 



1/2 



(3n — 2)mo 
2x~ 3 



1/2 



(A13) 



where m = mo is the reduced enclosed mass m(x) = 
ax 2 (nx — v) as x — > + , for a point mass at the centre. 

By free-fall solution (1A13[) and transformation ()A6|I 
with n + 7 = 2, we have a variable central mass accretion 
rate M given by 

M = k 3/2 t (3 - 3 ~ t) m /G . (A14) 

Only for the isothermal case of 7 = 1 (n = 1), can we have a 
constant M. We solve nonlinear ODEs (|A9[) and (|A10[) nu- 
merically using the fourth-order Runge-Kutta method. The 
SPS solution ()A11|I and asymptotic conditions (|A12|I and 
(|A13[) are all used to determine self-similar dynamic flows. 
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Figure Bl. A schematic illustration of the inside-out collapse 
scenario for the formation of double-peaked molecular line pro- 
files. The cloud consists of an outer static envelope and an inner 
collapsing core, separated by a spherical surface at r = r; n f a n 
in radial expansion. The loci (i.e. the two connected ovals) of a 
particularly chosen constant projected line-of-sight velocity com- 
ponent are shown. Two lines of sight are drawn towards right, 
labelled by 'Line of Sight 1' and 'Line of Sight 2', respectively. 
Line of Sight 1 passes through the core of central infall as well as 
the outer static envelope, that intercepts the ovals at points Bi 
and B2 for the same approaching projected velocity components 
and at points Ri and R2 for the same receding projected velocity 
components. Sufficiently far away from the centre, Line of Sight 
2 passes through only the outer static envelope. 



We also pay attention to the sonic critical curve (see Lou & 
Gao 2006, Lou & Wang 2006 and Hu & Lou 2008 for details) 



APPENDIX B: FORMATION OF 
DOUBLE-PEAK MOLECULAR LINE PROFILES 

Extensive spectral observations often show two peaks in high 
excitation emission molecular line profiles from star-forming 
molecular globules, and the blue peak is usually stronger 
than the red peak (i.e. blue asymmetric lines). This phe- 
nomenon can be explained in the framework of the inside- 



out collapse model i|Shulll977l ). We recount bel ow the rele- 
vant molecular line profile formation processes l|Zhou et al.l 
1 19931 ). 

Consider a spherically symmetric and non-rotating 
m olecular clo ud. In the isothermal inside-out collapse model 
oflShull| 19771) , the cloud consists of two radial regions, viz. (i) 
an inner collapsing core and (ii) an outer static envelope, as 
depicted in Fig. IB1I For simplicity, we use a static envelope 
here. The dynamical infall at the inner region determines the 
velocity profiles and hence the families of loci with equal 
magnitudes of projected line-of-sight velocity components. 
In Fig. IB1I we show two loci of the families that have equal 
magnitudes of projected line-of-sight velocity components 
but opposite directions. In the limit of no absorption (i.e. 
the optically thin regime), molecular lines from the emitters 
on these two loci will split, forming two symmetric com- 
ponents of peaks with equal but opposite (blue and red) 
frequency shifts. 

Let us first pick up Line of Sight 2 which passes through 
only the outer static envelope. This is a simple situation that 
the line has no bulk velocity broadening and splitting, i.e. a 



single-peaked line profile characterized by thermal and mi- 
croscopic turbulence broadenings would be expected. Next 
we consider Line of Sight 1 which passes through both the 
outer static envelope and the infalling core. Line of Sight 1 
intercepts the two loci of each family of constant projected 
velocity component at four points, namely Bi, B2, Ri and 
R2, with B and R implicating blue and red shifts respec- 
tively. We may denote the intensities of emission line compo- 
nents at these blue and red shifts as /blue and J re d- Suppose 
that the absorption of the emission line by the continuum 
is insignificant in the collapsing core. For illustration here, 
we may ignore the extinction in the envelope without los- 
ing generality, as it introduces only a scaling factor equally 
applied to both /blue and J rc d. With this simplification, the 
emission intensities are given by two transfer equations, 



S*i(l 



(Bl) 



/red = S 2 (l - e~ T2 )e- T1 + Si(l ~ e~ T1 ), (B2) 

(e.g. lZhou et al.lll993h . where in the vicinity of points 1 and 
2, ri and T2 are the optical depths, while Si and 5*2 are the 
source functions, with subscripts 1 and 2 denoting the outer 
and inner portions within the collapsing core at which four 
points (Bi, Ri) and (B2, R2) are located. It follows that the 
difference between the blue and red components of the line 
at this projected velocity component is 



/blue 



(Sa-ftXl-e-^Xl-e-^) 



(B3) 



For a local thermal equilibrium (LTE), the source func- 
tion S is given by th e local Planck function (e.g. 



iRvbicki fc Lightmanlll979l ). In the Rayleigh- Jeans regime of 
lower frequencies, the spectral intensity I v may be expressed 
in terms of the brightness temperature Tb in the form of 
Iv = 2u 2 ksTb/c 2 , where u, c and ks are the frequency, the 
speed of light and Boltz mann's constant, respectively (e.g. 
iDvson fc Williamslll997r i. Hence the difference in brightness 
temperatures between the blue and red peak emissions of a 
molecular line profile is simply 



Tb, rcd = (T a -Ti)(l-e~ T2 )(l 



(B4) 



Equation (IB4|I shows that a maximum contrast in the two 
peaks happens in the regime of n > 1 and T2 3> 1, 
and this gives a brightness temperature peak difference 
Tbiuc — Trod ~ T2 — T\ . Asymmetric line profiles arise only 
when there are substantial line optical depths in the hemi- 
spheres of the collapsing core that give rise to blue and red 
shifts of the line. In the opposite regime of ri <C 1 and 
T2 -C 1, the brightness temperature peak difference becomes 
Tbiuc — Trcd ~ (T2 — 2i)tiT2 ~ 0. In this case, the equal 
strengths of the two peaks can be understood as a direct 
consequence that in the absence of absorption (i.e. an op- 
tically thin gas) the emission from any location along the 
line-of-sight is visible. 

Molecular emission lines commonly used for cloud 
structure diagno stics (e.g. molecules H2CO, HCO + , CS, 

N2H 4 " etc. see Zhou. E vans & B utner 1990l ; Zhou et al 
— - — ■ 1 ^ — 



19931 ; IDvson fc Williamdll997l;IWard- Thompson fe Buckley 
200 it iLee. Bergin fc Evans! 12004 iTsamis et all 120081 ) are 
opaque in molecular globules that form protostars. For a 
collapsing core under the self-gravity, an increase in the tem- 
perature towards the centre, i.e. T2 > Ti, would then lead 
to a stronger blue peak than the red peak in the emission 
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Figure CI. An illustration of the computational grid in our ray- 
tracing radiative transfer calculations (with 4 3 = 64 cubic cells). 
The intensity of a molecular line is calculated cell by cell sequen- 
tially along each ray parallel to the z axis. The spectral image is 
formed in the (x, y) plane of sky. In our calculations for molecular 
emission line profiles, typically 101 3 cubic cells are used. 



line profile. The observed stronger blue peaks in the H2CO 
and CS line profiles from the molecular globule B335 thus 
support the notion of a collapsing core and in particular, 
the inside-out collapse scenario. But such clouds cannot be 
isothermal. Moreover, the temperature should be higher to- 
wards the cloud centre. 



APPENDIX C: RADIATIVE TRANSFER 
CALCULATIONS OF MOLECULAR LINES 

The ray-tracing radiative transfer equation in the LTE ap- 
proximation reads 



-(k c + K\)pl v + K c pS c + KlpS\ 



(CI) 



ds 

(e.g. IChandrasekhari 1960) where ds is the ray path differen- 
tial element, p is the mass density, S c and Si are the source 
functions, and k c and k\ are the mass absorption coefficients 
(or opacities), and the subscripts 'c' and T refer to associ- 
ations with the continuum and line emissions respectively. 
The corresponding optical depths are defined by 



ds' pn c (s') 



and 



Tl = 



ds'pmis') . (C2) 



For spectral line emissions, the opacity takes the form of 
where ki is the line centre opacity and dimen- 
sionless <j>„ is the normalized spectral line-profile function in 
radiation frequency v. The continuum remains more or less 
constant in the narrow frequency range of a spectral line 
profile; thus a molecular line profile is mainly determined 
by the opacity, source function and the intrinsic line profile 
function <f> u . For simplicity, we ignore continuum contribu- 
tions by setting k c — and S c — at this stage. Hereafter 
and in the main text, subscript T is also suppressed. 

To solve radiative transfer equation (IC1|I . we apply a 
forward ray-tracing computational algorithm on a three- 
dimensional (3D) Cartesian (x,y,z) grid with N 3 cells with 
N being the number of grids along each side. A spherical 
molecular cloud is embedded inside such a cubic grid. A ray 
propagates parallel to the z— direction towards an observer, 



and the spectral image of a cloud is projected onto the (x, y) 
plane of sky perpendicular to the z— direction (see Fig. IC1|) . 
We use N 3 = 101 3 cubic cells with an equal edge length 
of As for each grid cell which provide a sufficient spatial 
resolution for the purpose of our line profile computations. 

Approximately, each cell is locally homogeneous in den- 
sity, temperature and velocity, with their values taking the 
local values of a model molecular cloud at the cell centre. 
The numerical difference scheme for the local solution to 
transfer equation (|C1|I in a cell can be cast as 

- AT °+5 n (l-e^ AT ") , (C3) 



In = In-ie 



where I n is the emerging emission from the n-th cell, and 
5 n and Ar n are the source function and the optical depth 
of the n-th cell. By equation (|C2j) . the optical depth of the 
n-th cell is 



Ar n = pnAs = pKo(p(v, v p )As 



(C4) 



where p is the mass density of the cell, k is the effective 
line opacity, ko is the line-centre opacity and <j>(v, v p ) is the 
intrinsic normalized line profile. <j>(v, v p ) in equation (|C4[I 



is a function of velocity components v and v p , which cor- 
respond to the emission frequency v and the gas velocity 
in the cell projected along the line-of-sight respectively. We 
assume for simplicity that molecular line profiles are intrin- 
sically Gaussian with a normalized line profile function in 
the form of 



(2ttAu 



-1/2 



exp 



2Av 2 



(C5) 



where v p is the velocity component of the emitter projected 
along the line-of-sight and Av 2 is the square of the 'intrin- 
sic' characteristic broadening of a tracer molecular line in the 
absence of large-scale systematic dynamic flows. This intrin- 
sic broadening includes both the thermal and microturbu- 
lence effects. The source function of the n-the cell is given by 
the local Planck function (2/u//c 2 ){exp[/u//(fc s T)] - 
where h is the Planck constant. For emission lines in the 
Rayleigh- Jeans regime of low frequencies, we then have 
S n = 2v 2 kBT/c 2 1 where T is the thermal temperature of 
the cell. 

Our radiative transfer calculations are first carried out 
in a cell for a fixed frequency v and proceed sequentially 
to the adjacent cell along a ray parallel to the z— direction. 
Such calculations are then repeated for a range of frequen- 
cies along the same ray path. After finishing the integration 
along a ray, the same procedure will be repeated for all other 
rays on the grid, which then generates a molecular line spec- 
trum for each spatial pixel on the image plane. Depending 
on the desired spatial resolution, summation over the spec- 
tra of relevant pixels would give a line spectrum from a part 
of a molecular cloud. 



